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Abstract 

The hadronic light-by-light scattering contribution to muon g — 2 is examined 
using low energy effective theories of QCD, the Nambu-Jona-Lasinio model 
and hidden local chiral symmetry, as guides. Our result is —36 x 10~ n with 
an uncertainty of ±16 x 10~ n , which includes our best estimate of model de- 
pendence. This is within the expected measurement uncertainty of 40 x 10~ n 
in the forthcoming experiment at Brookhaven National Laboratory. Our re- 
sult removes one of the main theoretical obstacles in verifying the existence 
of the weak contribution to the muon g — 2. 
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The anomalous magnetic moment of the muon = | (g M — 2) is one of the basic physical 
observable that can be measured with high precision. The most recent results of experiment 
and theory are as follows 

a M (exp) = 1 165 923 (8.5) x KT 9 , 
a M (th) = 116 591 877(176) x KT 11 . (1) 

The uncertainty in the measurement of a M (exp) will be reduced in the new experiment at 
the Brookhaven National Laboratory to ~ 40 x 10~ n jj. This is about one-fifth of the 
one-loop weak correction |4j] 

a M (weak-l) = 195 (1) x KT 11 , (2) 

and of the same order of magnitude as the leading logarithmic term of the two-loop elec- 
troweak correction ||, a M (weak-2) = —42 x 10~ n , offering an exciting opportunity to test 
the quantum effect of the electroweak theory. 

Before comparing theory with the forthcoming measurement, however, it is necessary 
to reduce further the uncertainties in the theoretical prediction for the hadronic contribu- 
tion. The largest uncertainty comes from the hadronic vacuum-polarization contribution, 
a M (had.v.p.) |J. Fortunately, this contribution can be expressed as a convolution of known 
function with the experimentally measurable quantity R, the ratio of the hadron production 
cross section to the yU + /i~ production cross section in e + e~ collisions. Recent measurements 
of R at VEPP-2M will improve this estimate significantly @. Together with future mea- 
surements at DA$NE, BEPC, etc., the error in a M (had.v.p.) will be reduced to the level of 
the upcoming experimental limit. 

The contribution of the hadronic light-by- light scattering diagram shown in Fig. 1 is 
potentially a source of more serious difficulty because it cannot be expressed in terms of 
experimentally accessible observables and hence must be evaluated by purely theoretical 
consideration. Recently, some doubts have been raised 0,[| about the reliability of previous 
estimate [nj. In view of its importance in interpreting the experiment and in drawing 
inferences about potential "new physics" , we have decided to reexamine its theoretical basis. 



The bulk of hadronic light-by-light scattering contribution to a M is determined by the 
dynamics around the muon mass m M ~ 106 MeV, which is right in the midst of the non- 
perturbative regime of QCD. What we need is a reliable evaluation of an off-shell four-point 
function at these energies. In view of the current status of QCD, however, it is not an easy 
job to carry out such a calculation from first principles. 

Fortunately, this energy region is populated mostly by pions, and a considerable informa- 
tion is available about low energy pion dynamics. Chiral symmetry governs most of it. How- 
ever, higher energy regions may also contribute significantly to a M . Momentum-expanded 
amplitudes obtained in a systematic chiral expansion, however, cannot be introduced di- 
rectly into Feynman graphs for muon anomaly since it leads to divergent integrals. To get 
around this difficulty, Ref. used the Nambu-Jona-Lasinio (NJL) model. We chose to 



rely on the hidden local symmetry (HLS) approach fL2 |. Since the HLS Lagrangian can be 



derived from the extended NJL ( ENJL ) model fL3| , these two approaches are equivalent 
as far as application to our problem is concerned, 
(a) Relevant diagrams: 

The NJL model and the l/N c expansion suggest three distinct contributions to the light- 
by-light scattering amplitude at low energies JTJJ. Their contributions to are shown in 



Fig. 2. Fig. 2(a) shows the charged pseudoscalar meson loop contribution. It is 0(1) in 
1/N C expansion and 0(p 4 ) in chiral perturbation. Fig. 2(b) shows one pion pole diagram 
which is 0(N C ) and 0(p 6 ). Fig. 2(c) shows a quark loop diagram which is 0(N C ) and starts 
at 0(p s ) in chiral expansion. 

From the viewpoint of QCD, the single-quark-loop diagram Fig, 2(c) represents the 
averaged hadronic continuum effect in a certain energy region. At low energies, there should 
be higher-order QCD corrections, which can be approximated by a pion loop diagram Fig. 
2(a) and the Nambu-Goldstone boson pole diagram Fig. 2(b). Note that the latter two 
exclude a single-quark-loop contribution since the pion loop requires at least two quark loops 
and the pion pole starts from a diagram in which at least two gluons propagate between a 
quark and an antiquark forming the pion. 



Contributions involving more loops will be suppressed by a factor ?n At /(47r/ 7r ) ~ 1/11. 
Diagrams in which the pion loop in Fig. 2(a), tt° pole in Fig. 2(b), and u, d quarks in Fig. 
2(c) are replaced by kaon loop, 77 pole, and s quark, respectively, may also be nonnegligible. 

(b) Imbedding the light-by- light scattering amplitude: 

When the light-by-light scattering amplitudes are included in the Feynman graphs for 
g — 2, photons must be taken off-shell. The validity of low energy approximation to these 
amplitudes will then be affected by the convergence of photon momentum integration. 

In order to see what kind of problems might be in store, let us begin by comparing two 
previous treatments of the contribution to due to the charged pion loop diagram of Fig. 
2(a) ||10|| . One deals with the point-like 77T7T and 777T7T couplings (namely the scalar QED), 
which yields 



a N 3 



7T 



a^(sQED) = -0.035 57 (18) - . (3) 



(This is a slight improvement in statistics over that of Ref. fllP]). The second approach 



attempts to improve (H) based on the vector meson dominance model in which photons 
couple to charged pions through p meson. In Ref. [ID], this was achieved by modifying the 
photon propagator as 

i i m 2 p _ i i 



This leads to 



q 2 q 2 m 2 — q 2 q 2 q 2 — m 2 
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a M (nVMD) = -0.01 25(19) - . (5) 



Comparing (0) and @ we see that the introduction of the p meson makes a big differ- 
ence. This looks alarming. After all, vector mesons hardly contribute to tttt scattering near 
threshold since the chiral symmetry demands that their contribution vanishes in the soft 
pion limit. Why should they make such a difference, then ? Actually, this question is not 
posed correctly since the role of the p meson in (|5]) is primarily that of modifying the photon 
propagator and not directly related to the it-it interaction. What is more important is to 
check whether the vector meson contribution is properly included in (R). 



It was pointed out in Ref. || that the prescription given by (f|) does not respect the 
Ward identities for the coupling of photons to charged pions, even though it maintains 
gauge invariance. Upon closer examination, we found further that the substitution (^) is 
not consistent with chiral symmetry. 

In this paper, we solved these problems based on the HLS approach |T2|]. This is a 



convenient way to introduce the dynamics of pions and vector mesons preserving chiral 
symmetry and gauge invariance. It should be noted that this approach reproduces all current 
algebra results, such as the KFSR relation, as low energy theorems. In this sense, it is the 
leading candidate for the extension of chiral dynamics of pions to include vector mesons. 
Actually, insofar as only the low energy dynamics is relevant for the computation of a^, any 
model consistent with chiral symmetry should yield a similar contribution to a^. 

Before presenting our result, some comments relevant for computing are in order. 

1. The most important feature of the HLS Lagrangian, in computing a M (a) which corre- 
sponds to Fig. 2(a), is that it does not have the p°p 7r + 7r~ coupling. The naive substitution 
(|j) assumes the presence of this coupling in the Lagrangian. Indeed we found this to be the 
source of the problem pointed out in ||. Ward-Takahashi identity is satisfied once this is 
corrected. 

2. In computing 6^(6), the naive VMD model adopted in Ref. (JToJ did not have strong 



theoretical basis beyond that it provided an effective UV cut-off. It has been shown recently, 
however, that it is justifiable within the HLS approach f[4 |, at least as far as going off-shell 



with respect to photon momenta is concerned. This is also realized in the ENJL model, in 
which the ir° pole diagram contains two triangle loops of constituent quarks and p meson is 
allowed to propagate before the quark couples to the photon. 

3. To compute a^c), it is necessary to know how the quark couples to the photon and 
how vector mesons come into the picture. In this respect we are guided again by the ENJL 
model in which a quark loop couples to a photon through a vector meson (see Fig. 4 of Ref. 

II)- 

(c) Numerical results: 



An extensive numerical evaluation of the contribution of Fig. 2 to a M , within the frame- 
work of HLS, has yielded 

a M (o) = -0.003 55 (12) (^) 3 , 



a Jb) = -0.026 94(5) (- 

\7T 



3 
3 



a M (c) = 0.007 72(31) . (6) 

The errors quoted above are those of numerical integration only and do not include esti- 
mates of model dependence. The result a^a) is much smaller than (H). Similarly, a M (c) is 
considerably smaller than the corresponding results without vector mesons. a M (6) has a sign 
opposite to that of Ref. []10| |, which had a sign error in some terms. Thus, the difference 



between (|) and fl5|), which worried us a great deal, does not go away in spite of the im- 
proved theory which preserves chiral symmetry and Ward-Takahashi identity. This forced 
us to examine which regions of momentum space dominate in a M . We explored this problem 
by varying masses instead of examining momentum dependence directly, 
(d) Mass dependence: 

The dependence of a M on the mass m n or m q of internal loop in the light-by-light scat- 
tering amplitude were found to be as follows: 



aJa)(xm n ,M p ) ~ 2.81 x 10 2 x x 1 - 

,7T 



a 
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a M (6)(OTv,M p ) ~ -9.57 x 10~ 2 x x~ 2 (^j , 

a M (c) (xm q , M p ) ~ 1.14 x 10" 4 x x" 3 ' 7 3 , (7) 

for x > 3, where x is a scale factor of pion mass in a M (a) and a^b), or of quark masses (0.3, 
0.3, 0.5, and 1.5 GeV for u,d,s, and c, respectively) in a M (c). The first result in (|7D shows 
that the contribution of pion loop momenta drops off as x~ 2 as x increases. For instance, 
the contribution of pion momenta higher than 800 MeV accounts for only 7 percent of (0). 
From these results, we conclude that the hadronic light-by-light scattering amplitude, even 
when it is inserted in Feynman graphs for muon anomaly, can be described reasonably well 
by the graphs we have studied. 



We find the behavior a M (c) ~ x~ 4 quite encouraging. The fact that Fig. 2(c) contributes 
at all, as is seen from (|6|), implies that energy scale of (9(1 GeV) is important. On the other 
hand, the steep x dependence of a M (c) found in (|7|) is consistent with the fact that only the 
physical degree of freedom (mainly pions) is important at low energies 

We have also studied the dependence of a M on the vector meson mass: 

~m) 

a fl (b)(m n ,M) ~ [-0.069 3 + 0.31 



A* 

v7T. 

\ , / ( \ N 3 

M 



o„(c)K,M) ~ [+0.044 - 0.43 (^J] (-J (8) 

for M > M p . This shows that these functions decrease very slowly for large M. Such an M _1 
(instead of M~ 2 ) behavior seems to cast some doubt on the effectiveness of our approach 
since it implies that an appreciable contribution to a M comes from photon momenta far off 
mass-shell. Actually, this is entirely compatible with the dominance of low energy states 
in the light-by- light scattering amplitude. The naively expected M~ 2 behavior cannot be 
justified unless the M p — > oo limit and subintegrations in the Feynman diagram commute. 
As it is clear from power counting, this is not the case. Also, the results (^) show that there 
are strong cancellations between the first and second terms. In particular the cancellation 
in a M (a) for M = M p is almost complete, 
(e) Error estimates 

As stated above, three diagrams shown in Fig. 2 seem indeed to dominate the light-by- 
light scattering amplitude in the Feynman graphs for the muon anomaly. The contributions 
of kaon loop to a M (a), r] pole to a M (6), and heavier quarks to a^(c) can be readily included. 
We expect the error coming from further additional diagrams, as well as from the double 
counting possibility, to be less than the error caused by the approximate treatment of photon 
propagators in the HLS approach. By using HLS, we have not taken account, for example, 
of the continuum states above the vector mesons [|16|]. As long as we can restrict ourselves to 
pseudoscalar and quark loop diagrams, however, it is hard to imagine that these continuum 
states are relevant. Indeed the pion form factor is saturated with p meson. 



Taking these considerations into account, we estimate that the model-dependent errors 
from the terms of (|6|) are well within 20% of the M p -dependent second term. Including the 
contributions of the kaon loop coupled to in a M (a), the 77 pole to a M (6), the strange and 
charm quarks to a^c), we are thus lead to 

3 



a M ( a ) = -0.003 6 (64) (^) 
a„(&) = -0.032 8(66) (~ 



a N 3 



7T 

a x 3 



7T 



o (c) = 0.007 7(88) - . (9) 



(f) Summary: 

Using chiral symmetry, 1/N C expansion, and hidden local symmetry as guides, we have 
found that the hadronic light-by-light scattering amplitude can be represented reasonably 
well as the sum of diagrams containing a charged pion loop, diagrams with a 7r° pole, 
and diagrams with a quark loop at energies below 1 GeV. Based on this observation we 
have computed the hadronic light-by-light scattering correction to the muon g — 2 due to 
three diagrams of Fig. 2. The integration over the photon momenta receives considerable 
contribution from the region where the photons are far off shell. Estimating that these high 
mass contributions should be well within 20% of the vector meson contribution, we have 
obtained the result @ which leads to 

(light-by-light) = -36(16) x 10~ n (10) 

This is within the error expected in the upcoming experiment. Based on our analysis, and in 
view of the progress in the measurement of R 0, we are quite hopeful that the next round 
of experiments will indeed verify the weak interaction correction to the muon anomaly. 
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FIGURE CAPTIONS 



Fig. 1. Hadronic light-by-light scattering contribution (shown by the shaded blob) to the 
muon anomaly. Solid line and dashed line represent muon and photon, respectively. 

Fig. 2. Diagrams which dominate the hadronic light-by-light effect on a M at low energies, 
(a) Charged pseudoscalar diagram in which the dotted line corresponds to 7r ± , etc. (b) One 
of the 7T° pole graphs, in which the dotted line corresponds to tt° and the blob represents 
the 7T77 vertex, (c) Quark loop contribution, where quark is denoted by bold line. 
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FIG. 2 



